function [ hxu,hyu,dvu,dxr,dxl,dyt,dyb ] = getSpacingU( i,j,Xu,Yu )

if (i==1)
    dxl = 0;
else
    dxl = Xu(i) - Xu(i-1);
end

if (i==size(Xu,1))
    dxr = 0;
else
    dxr = Xu(i+1) - Xu(i);
end
    
if (j==1)
    dyb = 0;
else
    dyb = Yu(j) - Yu(j-1);
end
    
if (j==size(Yu,1))
    dyt = 0;
else
    dyt = Yu(j+1) - Yu(j);
end


hxu = 0.5*(dxr+dxl);
hyu = 0.5*(dyt+dyb);
dvu = hxu*hyu;



end

